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ABSTRACT 

We present high-resolution rotation curves and mass models of 26 dwarf galaxies from ‘Local Irreg¬ 
ulars That Trace Luminosity Extremes, The Hi Nearby Galaxy Survey’ (LITTLE THINGS). LITTLE 
THINGS is a high-resolution (~6" angular; <2.6 kms“® velocity resolution) Very Large Array (VLA) 

Hi survey for nearby dwarf galaxies in the local volume within 11 Mpc. The high-resolution Hi observa¬ 
tions enable us to derive reliable rotation curves of the sample galaxies in a homogeneous and consistent 
manner. The rotation curves are then combined with Spitzer archival 3.6/rm and ancillary optical U, 

R, and Uimages to construct mass models of the galaxies. This high quality multi-wavelength dataset 
significantly reduces observational uncertainties and thus allows us to examine the mass distribution 
in the galaxies in detail. We decompose the rotation curves in terms of the dynamical contributions 
by baryons and dark matter halos, and compare the latter with those of dwarf galaxies from The Hi 
Nearby Galaxy Survey (THINGS) as well as AGDM Smoothed Particle Hydrodynamic (SPH) simu¬ 
lations in which the effect of baryonic feedback processes is included. Being generally consistent with 
THINGS and simulated dwarf galaxies, most of the LITTLE THINGS sample galaxies show a linear 
increase of the rotation curve in their inner regions, which gives shallower logarithmic inner slopes a 
of their dark matter density profiles. The mean value of the slopes of the 26 LITTLE THINGS dwarf 
galaxies is a = —0.32 ± 0.24 which is in accordance with the previous results found for low surface 
brightness galaxies [a = —0.2 ± 0.2) as well as the seven THINGS dwarf galaxies (a = —0.29 ± 0.07). 
However, this significantly deviates from the cusp-like dark matter distribution predicted by dark- 
matter-only AGDM simulations. Instead our results are more in line with the shallower slopes found 
in the AGDM SPH simulations of dwarf galaxies in which the effect of baryonic feedback processes is 
included. In addition, we discuss the central dark matter distribution of DDO 210 whose stellar mass 
is relatively low in our sample to examine the scenario of inefficient supernova feedback in low mass 
dwarf galaxies predicted from recent AGDM SPH simulations of dwarf galaxies where central cusps 
still remain. 

Subject headings: Galaxies: dark matter halos - galaxies: cosmological N-body-l-SPH simulations - 
galaxies: kinematics and dynamics 


Electronic address: se-heon.oh@icrar.org 

^ The International Centre for Radio Astronomy Research 
(ICRAR), The University of Western Australia, 35 Stirling High¬ 
way, Crawley, Western Australia 6009, Australia 

^ ARC Centre of Excellence for All-sky Astrophysics (CAAS- 
TRO) 

^ Lowell Observatory, 1400 West Mars Hill Road, Flagstaff, Ari¬ 
zona 86001 USA 

^ Centre for Astrophysics Research, University of Hertfordshire, 
College Lane, Hatfield, ALIO 9AB United Kingdom 

^ IBM T. J. Watson Research Center, PO Box 218, Yorktown 
Heights, New York 10598 USA 

® Max-Planck-Institut fiir extraterrestrische Physik, Giessen- 
bachstrasse 1, 85748 Garching, Germany 

^ Cahill Center for Astronomy and Astrophysics, California In¬ 
stitute of Technology, 1200 East California Blvd, Pasadena, CA 
91125, USA 

^ Max-Planck-Institut fiir Astronomic, Konigstuhl 17, 69117 
Heidelberg Germany 

^ National Radio Astronomy Observatory, 1003 Lopezville Road, 
Socorro, NM 87801 USA 

Physics Department, New Mexico Institute of Mining and 
Technology, Socorro, New Mexico 87801 USA 

Department of Physics, Florida International University, CP 
204, 11200 SW 8th St, Miami, Florida 33199 USA 

Australia Telescope National Facility, CSIRO Astronomy and 
Space Science, PO Box 76, Epping NSW 1710, Australia 

Department of Physics, Pennsylvania State University Mont 
Alto, Science Technology Building, Mont Alto, PA 17237 USA 

School of Physics and Astronomy, University of Southampton, 
Southampton S017 IBJ UK 


1. INTRODUCTION 

Dark matter the existence of which is indirectly in¬ 
voked by its gravitational effect in individual galaxies 


dark energy, the energy budget in the Universe (Zwicky 

1937 

van den Bergh||1961[ |Rubin & Ford| 

19701 [Bosma 

1978 

Peebles||1982l IRiess et al. 119981 Per 

mutter et al. 

1999 

Spergel et al. 2003). In particular, the role of 


evolving galaxies in the early Universe but also in shap¬ 
ing the lar ge-scale structure in the Universe through cos- 
mic time (Blumenthal et al.l 1984 Colless et al. 2001a 


Padmanabhan et aL| |2007[ [Jones et aL| |2009[ |Komatsu| 


et al. 


2011 etc.). The cosmological importance of DM 


has driven efforts to explore the physical nature of DM 


particles and to attempt their direct detection ([Moore 


et al.|[2Q0T| Gaitskell| 2004 Angloher et al.|[2012[ [Akerib 

et akf 20141 and references therein). Of many candi- 
dates for DM particles, Cold Dark Matter (CDM) has 
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been envisaged as one of the most successful models in 
that numerical simulations based on a paradigm combin¬ 
ing CDM with the cosmological constant, A (so-called 
ACDM). This describes well the large-scale struct ure in 


(Ott 




et al.l 12012 

TTKIGS^l 


2012 

), LITTLE THINGS22 ( 

Hunter et al. 


Begum et al. 

2008 

LSHJ 

[EEB^l 

Gan- 

2011) LVHlb^lKdnra 

ski 

2010 

), have al- 


the Universe traced by surveys such as the SDSS^^ ( York 
Doroshkevich et al. 2004), 2dFGRS^^ ((Jol- 


et al. 2000 Dor oshkevich e t al. 2004), 2dFG RS^° ((Jol- 
less et al.|200lb|), bdhOy^y dJones et al.|2004l), and CMF 
power spectrum observations dPrimack ( 2003 _ Spergel 


et al. 2003, 2007 Planck Collaboration et al. 2014). 


However, despite the success ot cosmological ACDM 

simulations in producing the large-scale structure of the 
Universe, distinct differences between the simulations 
and observations have been found in the DM distribu¬ 
tion near the centre of individual galaxies. The simu¬ 
lations have consistently predicted cusp-like DM distri¬ 
butions towards the centers of galaxies, being described 
by a pow er law, p^R°" where R is the galaxy radius and 


non et al. 

lowed us to derive more reliable Hi rotation curves of 
galaxies and examine their central mass distributions 
within 1 kpc where the predictions of ACDM simula¬ 
tions are most distinctive. For example, high-resolution 
(0.1-0.2 kpc) DM density profiles of seven dwarf galaxies 
from THINGS, complemented wi th the ^Spitzer Infrare d 
Nearby Galaxie s Survey’ (SINGS ;|Kenn icutt et al.||2003|), 
were derived by Oh et al.l (2008 |2011b ). hfom this, they 
found that the mean value of the inner density slopes, 
a, of the seven dwarf galaxies is —0.29 ± 0.07 which is 
in good agreement with the value of —0.2 ± 0.2 derived 
earlier from a la rger number of Low Surf ace Brightness 


a~—1.0 (Moore 1994[ Navarro et al. 


1996b 


et al.| |1997 ' Moore et al.| |1999[ [Ghigna et a 


Navarro 


Klypin et al 

120011 IPower et al.l 12003| IStoehr et al. 

20031 Navarro et al.| 2004a 

diemand et al.l 12008 Stadel 

et al.| 120091 INavarro et al. 

20101 llshiyama et al. 2013 


etc.). In contrast, inferred DM distributions in nearby 
dwarf galaxies in the local Universe have shown a linear 
velocity increase towards their centers, gi ving rise to a 
sizable density-core io 


de B 

ok et al.ll 

1996[ |de Blok & McGaugh||19y7| |de Blok| 

et al. 

20011 de Blok & Bosma|20021 1 Weldrake et al. 

2003 

Spek 

tens et al. 20051 Kuzio de Naray et al. 112006 | 

<(uzio 

de Naray et al.l 120081 lOh et al.l 120081 lOh et al.l 2011bl 


(LSB) galaxies (de Blok & Bosma 2002). In the past, 
it has been argued that observational systematic effects 
such as beam smearing, center offsets and non-circular 
motions could have affected the derived central DM dis¬ 
tributions of galaxies as these observational biases tend 
to flatten the derived inner DM density profiles, hiding 
the central cusps. However, because of the high reso¬ 
lution and quality of the above mentioned observations, 
these potential biases were significantly reduced. As a re¬ 
sult, observational evidence for the core-like distribution 
of DM near the centers of dirr galaxies is particularly 
strong. 

In order to explain this behaviour, baryonic feedback 
processes have been proposed as a means for removing 
the central cusps expected from ACDM DM-only simu¬ 
lations. More specifically, it is expected that DM and 


a discussion on gas kinematics which produces shallower 
density profiles than those from stellar kinematics). This 
clear discrepancy of the central DM distribution in galax¬ 
ies between AGDM simulations and observations, the so- 
called ‘cusp/core’ problem has been one aspect of the 
small-scale crisis in ACDM cosmology whic h is likely con¬ 


necte d to the ‘missing satellites’ problem (Brooks et al. 


20131 and ‘too-big-to-fail’ proble m (dense satellites; see 


frequent explosions of supernovae (SNe) (Larson 

1974 

Navarro et al.lll996allDekel et al. 1120031 IMo & Mao 

MM 

Mashchenko et al.l 20081 de Souza et al. 2011 . 

3rook 

et al. 

20111 IPontzen & Governatoll2012| IDi Gintio et al. 

2014 

Oglya k Morl||20l4 

etc.). However, due to numer- 


Boylan-Kolchin et al. 2011 


2012). Given that DM com- 


bined with the AGDM paradigm is successful in explain- 


well as the lack of understanding of the detailed baryonic 
physics in galaxies, taking the baryonic feedback into 
account in hydrodynamical simulations of dwarf galax¬ 
ies was considered difficult and limited to simulations of 


as galaxy formation and evolution, there are good rea- 

high-redshift galaxies (e.g., Mashchenko et al. 

2006) 

sons to explore ways to resolve the ‘cusp/core’ problem. 

Governato et al. (2010 

) were the hrst to peri 

orm h 


(2014) for the latest review of the ‘cusp/core’ problem. 


Late-type dwarf galaxies in the local Universe, with 
a simple dynamical structure (no bulge and spiral com¬ 
ponents) have been used for addressing the central DM 
distribution in galaxies since their dynamics are usually 
dominated by DM, enabling us to derive more accu¬ 
rately the DM distribution near their centers. These 
diffuse dwarf galaxies have provided significant obser¬ 
vational constraints on the central DM distribution in 
galaxies. Over the past decade, several high-resolution 
neutral hydrogen (Hi) surveys of galaxies in the local 
Universe (<11 Mp c) using radio interf erometers, among 


others THINGS^o (Walter et al. 2008), VLA-ANGST^i 


The Sloan Digital Sky Survey (http://www.sdss.org/) 

The 2dF Galaxy Redshift Survey (http://www2.aao.gov.au/ 
2dfgrs/) 

The 6dF Galaxy Survey (http://oldweb.aao.gov.au/local/ 
www/6df/) 

The Hi Nearby Galaxy Survey (http://www.inpia-hd.inpg.de/ 


resolution cosmological N-body-|-Smoothed Particle Hy¬ 
drodynamic (SPH) simulations of dwarf galaxies which 
include the effect of detailed baryonic feedback processes, 
in particular physically motivated gas outflows driven by 
SN explosions. From this, they found that the photo¬ 
metric and kinematic properties of the simulated dwarf 
galaxies are in close, qualitative agreement with those of 
observed nea rby dwarf galaxies. More quantitatively, |Oh| 
et al. (2011a I performed an analysis of the baryonic and 
DM mass distributions of the simulated dwarf galaxies, 

THINGS/Dverview.html) 

Very Large Array - ACS Nearby Galaxy Survey Treasury 
(https://science.nrao.edu/science/surveys/vla-cuigst) 

Local Irregulars 'i'hat 'i'race Luminosity Extremes, The 
Hi Nearby Galaxy Survey (https://science.nrao.edu/science/ 
surveys/littlethings) 

h'aint Irregular Galaxies GMRT Survey 

The Survey of HI in Extremely Low-mass Dwarfs (http:// 
www.macalester. edu/'^jcannon/shield .html) 

The Local Volume Hi Survey (http://www.atnf.csiro.au/ 
research/LVHIS/) 







































































































































































































High-resolution mass models of dwarf galaxies from LITTLE THINGS 
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in exactly the same way as the THINGS dwarf galax¬ 
ies were analysed, and showed that their derived rota¬ 
tion curves and the corresponding DM density profiles 
are consistent with those of the THINGS dwarf galaxies 
and show a linear increase of velocity in the inner re¬ 
gion inherent of shallow DM density profiles. This sug¬ 
gests that repeated gas outflows driven by SN explosions 
even without a burst of star formation are able to play 
a fundamental role in removing the central cusps and 
inducing flatter DM density slopes near the centers of 
dwarf galaxies. This is in contrast to clusters of galaxies 
where galaxy interactions are more likely to be the dom- 


tral DM distribution of galaxies in Section Lastly, we 
summarise the main results of this paper ana conclusions 
in Section [71 

2. THE DATA 

We use high-resolution Hi data of 26 nearby (<11 
Mpc) dwarf galaxies from LITTLE THINGS to address 
the central DM distribution of the galaxies. The sam¬ 
ple galaxies show a regular rotation pattern in their 2- 
dimensional (2D) Hi velocity fields (see Appendix), al¬ 
lowing us to derive reliable rotation curves which include 
the contributions to their kinematics of both their DM 


than star formation activities (Richtler et al. 2011 

). See ing the distances of the sample galaxies, the linear res- 

Governato et al. 

2010) and Pontzen &: Governato 

2012 

olutions of the LITTLE THINGS Hi data (~6") range 
from ^26 to 200 pc with an average of 100 pc which is 
sufficient to resolve the inner 1 kpc region of the galaxies 
in so far unmatched detail. This enables us to examine 

2U14p lor detaifei 

explosions on the 
THINGS was c 

discussions regarding the effect of Si\ 
central cusps. 

inly able to probe a small number (i.e. 


seven) of dwarf galaxies. It is therefore essential to ex¬ 
tend the investigation to a larger number of dwarf galax¬ 
ies in order to obtain a statistically robust observational 
sample to which simulations can be compared to. Data 
on a larger sample of dirr galaxies has now been pro¬ 
vided by t he latest Hi survey of nearby galaxies, LITTLE 


THINGS (Hunter et al. 2012). LITTLE THINGS is a 
high-resolution (~6" angular; < 2.6 kms“^ velocity reso¬ 
lution) Hi survey for 41 nearby (< 11 Mpc) gas-rich dwarf 
galaxies undertaken with the NRAO^® Very Large Ar¬ 
ray (VLA) in the northern sky. The Hi observations are 
complemented with other wavelen gth data, such as Ha, 
optic al [/, B, V, and near infrared (Hunter & Elmegreen 
20061, archival Spitzer infrared and GALEX ultraviolet 


images as well as follow-up observations with ALMA and 
Herschel. These high-quality multi-wavelength data sets 
significantly reduce the observational uncertainties inher¬ 
ent in low resolution data which may result in hiding the 
central cusps. Of the 41 galaxies, we select a sample of 26 
dwarf galaxies (three of them are also in THINGS) which 
show a regular rotation pattern in their velocity fields. 
In this paper, we extract (I) bulk and non-circular mo¬ 
tions of the sample galaxies, (2) derive rotation curves, 
(3) decompose the derived rotation curves in terms of the 
contributions by baryons and DM halos, and (4) address 
the central DM distribution by making a direct compari¬ 
son between the derived DM distributions of the galaxies 
to those of SPH-|-N-body simulations of dwarf galaxies. 

The structure of this paper is as follows. The data used 
for deriving the mass models of our sample galaxies are 
described in Section In Section ^ we present the ro¬ 
tation curves, the mass models of baryons, and the DM 
mass modeling of the galaxies. In Section [^ we discuss 
the central DM distributions of the sample galaxies by 
comparing them with those of dwarf galaxies from both 
THINGS and simulations. We then discuss the effect of 
SN feedback on the central cusp in Section followed by 
the discussion of the effect of beam smearing on the cen- 

NRAO is a facility of the National Science Foundation oper- 
ated under cooperative agreement by Associated Universities, Inc. 
These data were taken during the upgrade of the VLA to the Ex¬ 
panded VLA or EVLA. In this paper we refer to the instrument 
as the VLA, the retrofitted antennas as EVLA antennas, and non- 
retrofitted antennas as VLA antennas. This emphasizes the hybrid 
nature of the instrument and distinguishes it from the far more 
powerful Jansky VLA or JVLA it has become since 2012. 


the central DM distributions of the galaxies in detail. 
In addition, observational systematic effects inherent in 
low-resolution data (e.g., beam smearing, kinematic cen¬ 
ter offset and non-circular motions) are significantly re¬ 
duced in the high-resolution Hi data, which allows us to 
derive more accurate underlying kinematics of the sam¬ 
ple galaxies. 

Although the total kine matics of late-type dwarf galax- 


ies is dominated by DM (|Prada & Burkert 20021, it is 
nonetheless important to separate the contribution by 
baryons from the total rotation curve. This is achieved 
by using Spitzer archival IRAC O.Ojum and ancillary 


optical color information (Hunter & Elmegreen 2006). 
Spitzer IRAC 3.6/im images are much less affected by 
dust than maps at shorter wavelengths and trace the old 
stellar populations that occu py the dominant fract ion of 


the stellar mass in galaxies (Walter et al. 2007). Al¬ 


though there is some contarnin ation by pAH erni ssion 


extending into the 3.6p,m band (Meidt et al. [2014) this 
is reduced in dIrr galaxies which, by virtue of naving low 
heavy element abundances, have a correspondingly lower 
dust content. _ 

Acc ording to the simulations by [Governato et al.] 
(2012), the degree of baryonic feedback is largely depen- 
dent on the amount of stars, so past star formation ac¬ 
tivity in the sense that for the same total mass budget 
(baryons -I- DM), those galaxies with a higher past star 
formation activity will have had more significant out¬ 
flows. Therefore a reliable measurement of the stellar 
mass in a galaxy is essential to investigating the effect 
of baryonic feedback on the central cusp. The basic ob¬ 
servational propertie s of the sampl e galax ies are listed in 
Table [T| We refer to Hunter et al. (2012) for a complete 
description of the Hi observations and data reduction. 

3. MASS MODELS 

In this Section, we perform mass modeling of the 
26 LITTLE THINGS dwarf galaxies using the high- 
resolution VLA Hi data and Spitzer IRAC 3.6^m images 
as described in SectionAs mentioned earlier, the sim¬ 
ple dynamical structure and dominant circular rotation 
in the disk of the sample galaxies help to reduce un¬ 
certainties stemming from the decomposition of galaxy 
rotation curves into contributions due to baryonic and 
non-baryonic matter. 

The first step in the mass modelling includes deriv- 
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ing rotation curves of the sample galaxies using the 2D 
Hi velocity fields which reflect the total kinematics in¬ 
cluding both baryons and a DM halo. Hi is mostly dis¬ 
tributed in the disk of a galaxy where the circular ro¬ 
tation is dominant and is a useful kinematic tracer for 
deriving the galaxy rotation curve. This is mainly due 
to the larger radial extent of Hi in the disk compared 
to ste llar components (e.g., 3—4 times; Sofue & Rubin 
20011. We proceed to derive mass models of the baryons 
(gas and stars) using Hi integrated intensity maps and 
Spitzer IRAC 3.6/rm images, and subtract their contri¬ 
bution from the total kinematics. Lastly, we quantify the 
kinematic residuals in order to examine the DM distri¬ 
bution near the centers of the galaxies. In the following 
sections, we describe these mass modelling procedures in 
more detail. 

3.1. Rotation curves: total kinematics 
3.1.1. Tilted-ring fits 

For the derivation of the rotation curves of our sam¬ 
ple galaxies, we fit a 2D tilted-ring model which con¬ 
sists of a series of concentric ellipses to the 2D velocity 
fields ext racted from the Hi data cubes of the sample 
galaxies (Rogstad et al. |1974 ). Each ellipse has its own 
geometric and kinematic parameters, such as center po¬ 
sition (XPOS, YPOS), position angle (PA), inclination 
(INCL), systemic velocity (VSYS), and rotation velocity 
(VROT). This so-called ‘tilted-ring analysis’ of 2D ve¬ 
locity fields obtained from Hi, CO or Ha spectroscopic 
observations has been widely used for deriving rotation 
curves of disk-dominated galaxi es (e.g. , rotation curves 


of the THINGS galaxies sample; de Blok et al. 20081. 

Tilted-ring models, however, only hold for those cases 
where the velocity held is a reliable representation of 
the overall kinematics of a galaxy. Non-circular motions 
caused by star-forming activity like stellar winds, SNe, 
etc., as well as spiral arms, a bar-like or triaxial galaxy 
potential, galaxy mergers or tidal interactions disturb gas 
motions in galaxies on small and large scales, resulting 
in distorted velocity helds. 

Low resolution Hi data usually smooth any features re¬ 
lated to small-scale non-circular motions. This is known 
as beam smearing and tends to yield a velocity gradient 


ularl 

Y in the central regions (e.g., |de Blok & McGaugh 

1997 ISwaters et al.| 20001 

van den Bosch & Swaters|200Il 

McGaugh et al. 200ip. 

3eam smearing is signihcantly 
n data such as THINGS and 

reduced in high-resolutio 


LITTLE THINGS (see de Blok et al. 

2001 

et al. 

mi * 

2001 

* I im 

Kuzio de Naray et al. 

ill 1.^ . m/im - 1 — 

2006 

• 


In order to correct for small-scale random motions and 
extract only the component due to circular rotating ve¬ 
locity in a g alaxy, we derive the ‘bul k velocity held’ as 


proposed by Oh et al. (2008 2011b). Compared with 


other typical types of velocity helds, such as intensity- 
weighted mean (IWM), single Gaussian ht and hermite 
h3, the bulk velocity held has been found to be ideally 
suited to extracting the underlying bulk rotation of a 
galaxy in the presence of random non-circular motions. 
We point out, as an example, how the ‘kinks’ or ‘wiggles’ 
of the iso-velocity contours of the extracted bulk velocity 
helds displayed in panel (e) of Fig. A.l are weaker than 
those of the IWM velocity helds (panel b). The bulk 
velocity held appears to much better represent the over¬ 
all kinematics of the galaxy and hence the underlying 
gravitational potential. 

Fol lowing the standard procedure described in |Bege-| 


man 


(19891, we ht tilted-ring models to the bulk velocity 
helds of the samp le galaxies. For this, we u se the ‘rotcur’ 
task in GIPSY^^ (Ivan der Hulst et al.jl992 1. The derived 
rotation curves of the galaxies are presented in the hg- 
ures of the Appendix (e.g., A.2). As seen from the scatter 
in the hts made with all ring parameters free (open cir¬ 
cles) in the hgures, the extracted bulk velocity helds are 
not completely free from the effect of small-scale random 
motions in the galaxies. However, they are relatively in- 
signihcant and are averaged out after several iterations 
as shown in the hnal rotation curves (solid lines) in the 
hgures. 

This is also conhrmed in the harmonic analysis of 
the ve locity helds. As described in |Schoenmakers et al.| 
(19971, we perform harmonic decompositions of the bulk 
and IWM velocity helds of the sample galaxies. For this, 
we use the task ‘reswri’ in GIPSY. We expand the ve¬ 
locity helds into sine and cosine terms up to 3'^'^ order 
(i.e.. Cm and Sm where m = 1,2 and 3) after hxing the 
center position, PA and INGL with those derived from 
the tilted-ring analysis. If we allow RESWRI to ht a ve¬ 
locity held with center position, PA and INCL as free 
parameters, non-circular motions tend to be absorbed 
into variations in these geometrical parameters, underes¬ 
timating the amount of non-circular motions. As an ex¬ 
ample, streaming non-circular motions in a barred galaxy 
are mainly responsible for the radial moti ons that are 


in this study, with few exceptions, appear to show no 


typically rehected in the si and S 3 terms (Schoenmak- 
s et alT] |1997[ Wong et al. 2004[ Spek cens fc be 
wood||2007p. However, these radial motions can also be 
modelled by a radial variation of PA, without the need 
for Si and S 3 terms. Similarly, other ring parameters 
can affect the harmonic analysis in the same way if they 
are kept as free parameters in the ht. To quantify and 
describe non-circular motions, we calculate the absolute 
amplitudes (A) and the phases of each component de¬ 


signihcant large-scale kinematic features, such as bars, composed. Eor the am 

soiral arms or warns in their Hi or Svitzer IRAC 3.6/im Am.(R) as described in 

plitudes, we take the median of 

Schoenmakers et al. 

(1997) (see 

images. However, our sample galaxies do suffer from the also 

Trachternach et al. 



effect of small-scale turbulent gas motions due to stellar 
winds and SNe. In general, the derived velocity held in 
dwarf galaxies is more vulnerable to the impact o f stellar 
activity due to their low gravitational potential (Walter 
et al.| 1998[ |Walter fc Brinks||2001| see also! 


et al.||20Iip. Therefore, the extraction of a robust veloc¬ 


ity held is essential if one wants to derive the undisturbed 
underlying kinematics of a galaxy. 


Ai(i?) = ys^(i?P, 

for m > 1 , 

Am{R) = ^JCm{.Rf + Sm{RY, 

The Groningen Image Processing System 


( 1 ) 


(2) 
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where R is the galaxy radius. 

As shown in the section labelled ^Harmonic Analysis^ 
in the Appendix, the amplitudes of harmonic terms (e.g., 
C 2 , Si and S 2 which are corrected for inclination) derived 
from the bulk velocity fields (black dots in the Harmonic 
Analysis panels) are lower than those derived using the 
IWM velocity fields over all radii. This shows that the 
effect of random non-circular motions are largely reduced 
in the bulk velocity fields. The tilted-ring parameters 
of the sample galaxies derived using the bulk velocity 
fields are given in Table We note that the kinematic 
center positions given in Table[^that are derived from the 
tilted-ring analysis are offset fr om the optica l (H-b and, 
mostly) central isophot used by [Hunter et al. (2012|. 


3.1.2. Asymmetric drift correction 


Pressure support caused by random gas motions in the 
gaseous disk of a galaxy tends to lower the rotation veloc¬ 
ity, which results in an underestimate of the dynamical 
mass of the galaxy. In general the dynamical effect of the 
pressure support is higher in the outer region of a galaxy 
where the gas density is low. In particular, this, so-called 
‘asymmetric drift’ is significant in dwarf galaxies whose 
maximum rotation velocities are comparable to the ve¬ 
locity dispersions in the gas disk. The asymmetric drift 
correction hence should be made to derive more reliable 
rotation curves of galaxies where the dynamical support 
by random motions to its gas disk is significant. This is 
the case of some of our sample galaxies whose velocity 
dispersions are comparable to their maximum rotation 
velocities in the outer region. For the asymmetric drift 
correction , we fo llow the method describe d in Bureau ^ 
Carignan (2002 see also Oh et al. 2011b I as follows: 

'I'he asymmetric drift correction (Td is given as. 


: —Ra 


= —Ra' 


^ d\n[pa^) 

dR 

91n(EHio-^) 

dR 


(3) 


where ctd is the asymmetric drift correction, R is the 
galaxy radius, cr is the velocity dispersion, and p is the 
volume density of gas disk. In general, p can be ap¬ 
proximated as the gas surface density Ehi for a gas disk 
with an exponential distribution in the vertical distri¬ 
bution and a constant scale height zq (i.e., d(ln(z(^/dr 
= 0). Large fluctuations in the derivative in Eq. ^can 
be smoothed by fitting an analytical function whicm has 
three free parameters, /q [Mq pc“^ km^ s“^], i?o [arcsec], 
and a [arcsec“^] to the numerator as follows. 


Sma^(R) 


/o(i?o + l) 
Ro + e“« ■ 


(4) 


Lastly, the corrected rotation velocity 14or is derived by 
adding the asymmetric drift correction a^ to the rotation 
velocity Vrot derived from tilted-ring fits, quadratically, 

Kor = Kot+4- (5) 

The analytical function given in Eq provides a good 
fit to most sample galaxies except for DDO 52, IC 10, 
NGC 3738 and UGC 8508 where a significant degree of 
radial fluctuation is present in at small radii. How¬ 

ever, the correction at small radii is insignificant and thus 
will not affect significantly the final results. 


Valenzuela et al. (2007) have shown that the standard 


asymmetric correction can be underestimated without 
considering the gas pressure gradients triggered by star 
formation and feedback. However, the effect of additional 
thermal pressure gradients in the gas is most likely in¬ 
significant for the sample galaxies in this study. The 
thermal pressure is already included in the observed ve¬ 
locity dispersion used for the asymmetric drift correction 
as part of its broadening. In addition, as shown in the ve¬ 
locity dispersion map (i.e., moment 2) in the Appendix, 
significant anisotropy in the dispersion is not found in 
the galaxies. This implies that any separation of tur¬ 
bulent and thermal pressures, with explicit calculation 
of asymmetries in the turbulent component, would have 
insignificant effect on the asymmetric drift corrected ro¬ 
tation curves. The asymmetric drift corrected rotation 
curves of the sample galaxies which are used for the mass 
modelling are presented in the Appendix. 


3.2. Gas distribution 

The rotation curves derived in Section [3 .1 1 already pro¬ 
vide a good approximation of the DM halos’ kinematics 
of the sample galaxies given the dominant contribution 
of the DM halo to the total kinematics of dwarf galaxies. 
However, to derive more accurate DM distributions of 
the galaxies, we construct mass models of their gaseous 
and stellar components which account for most of the 
baryons in dwarf galaxies. 

We use total integrated Hi intensity maps (moment 0) 
of the galaxies to derive the mass model of the gaseous 
component. For consistency with the rotation curves in 
Section |3.H we apply the derived tilted-ring models to 
the Hi intensity maps, and obtain gas surface density 
profiles of the galaxies which are scaled up by a factor 
of 1.4 to take Helium and metals into account. We then 
convert the gas surface density profiles to the correspond¬ 
ing gas rotation velocities assuming that gas components 
are mainly distributed in a thin disk. As an example, the 
derived gas surface density profile and the corresponding 
rotation velocity of CVnIdwA are shown in the panels (g) 
and (h) of Fig. A.3, respectively. Here, we do not correct 
for the effect of molecular hydrogen (H 2 ) since low metal- 
licities in dwarf galaxies can induce only a small fr action 


of the gaseous component in the form of H 2 

(e.g., 

Leroy 

et al. 

2007 

Leroy et al. 

2011 

Schruba et al. 

2012 

1 . 


3.3. Stellar distribution 


We use Spitzer IRAC 3.6/im images to derive mass 
models of the stellar components of the galaxies. Gom- 
pared with optical images, the Spitzer IRAC 3.6/im im¬ 
age is less affected by dust and less sensitive to young 
stellar populations which usually emit most energy in 
the optical regime but occupy only a small fraction of 
the total stellar mass. Instead, the Spitzer IRAC 3.6/im 
image is useful for tracing old stellar populations that are 
dominant in late-type dwarf galaxies. This enables us to 
derive a robust estimate of the stellar mass of our sam¬ 
ple galaxies as used for deriving the mass models of the 
stellar components of THINGS galaxies dde Blok et a l. 
'20081 iTrachternach et al. 1 12008 lOh et al.l 20081 |2011bp. 


"Clke we did for the gas component, we derive 3.6/im 
surface brightness profiles of the stellar components of 
the galaxies by applying the derived tilted-ring parame¬ 
ters to the Spitzer IRAC 3.6/im images as shown in the 
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Fig. 1.— Comparison of the stellar masses derived usin g our method with those derived using the multi-band spectral energy distribution 
(SED) fitting technique described in|Zhang et al.| (|2012||. The solid and dashed lines indicate a least-squares fit with a slope of 0.96, and 
line of equality, respectively. The Icr scatter ot the tit to the data is 0.20 dex. 


figures in the Appendix (e.g .. Fig A.3). In general, 
as discussed in Walter et al. (2008), the Spitzer IRAC 
3.6/xin image provides a pseudo dust free picture of old 
stellar populations in galaxies. However, unlike the case 
of the gas component whose mass can be directly esti¬ 
mated from Hi observations, estimating the stellar mass 
in galaxies is critically dependent on the assumed stel¬ 
lar mass-to-light ratio (T*) value which usually gives rise 
to the largest uncertainty when converting the luminos¬ 
ity profile to the mass density profile. In order to de¬ 
rive more reliable T* values in the 3.6/im band, we use 
an empirical relation between galaxy optical colors and 
values based on stellar population synthes is models 
(Bruzual fc C harlot | 2003 Bell fc de Jong|[20M ) as given 
in Oh et al. "(20081. Hereafter, we call these as model 


Tj-® values. Using the derived ® values, we convert 
the 3.6/rm surface brightness profiles of our sample galax¬ 
ies to stellar surface density profiles (see the figures in the 
Appendix, e.g., Fig. A.3). 

As discussed in Oh et al. (2011a), the kinematic 


method combined with the model T values has been 
found to be reliable for estimating stellar masses of late- 
type dwarf galaxies based on a comparison of the derived 
stellar masses of the simulated dwarf galaxies to the input 
ones. This robustly supports the veracity of the method¬ 
ology used for measuring the stellar masses of our sample 
galaxies. In addition, we also refer to the stellar masses 
of the sample galaxies derived usin g a spectral energy 
distribution (SED) fitting technique (Zhang et al. 2012). 
As shown in Fig.[^ the stellar masses derived using these 
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two independent methods show good agreement within 
the scatter. 


Following Oh et al. (2011b), we calculate the corre¬ 
sponding rotation velocities of the stellar components 
of the sample galaxies from the derived surface density 
profiles. For this, we assume a vertical sech^( 2 ) scale 
height distribution of stars with a ratio of hjZQ=h where 
h and zg are the radial scale length and the vertical scale 
height of stellar diskTn_the_T6£tni£urfa£ebri£htaess_gro- 
files, respect ively (van der Kruit & Searle 1981| Kregel 


et al.| |2002|). Altnough it may overestirnate the rota¬ 
tion velocities of the stellar components of galaxies with 
a fatter stellar disk, this is a valid assumption for most 
disk-dominated dwarf galaxies like our sample galaxies. 
The derived rotation velocities of the stellar components 
of the sample galaxies are shown in the figures in the 
Appendix (e.g., panel (d) of Fig. A.3), and the stellar 
masses estimated in this paper are given in Table 1^ 

The mean ratio of the masses between the gas and stel¬ 
lar components of the sample galaxies, < Mgas/M^™ > 
is ~ 5.6 which is consistent with < Mgas/M®™ >~ 7.2 
(where Mf™ is derived using a SED fittin g technique) 
of 34 LITTLE THINGS sample galaxies in [Zhang et al. 
(2012). This indicates that the majority of the baryons in 


our sample galaxies is in the form of gaseous components. 
As discussed earlier, given that the mass of gaseous com¬ 
ponents in galaxies can be reliably measured from Hi 
observations without any critical assumption, the mass 
models of the baryons derived in this study are likely to 
provide a good description of the distribution of baryons 
in the galaxies. Therefore, any remaining uncertainties 
in the mass models of baryons are most likely to be in¬ 
significant and thus will not affect significantly the final 
mass models of the DM halos of our sample galaxies. 


3.4. Disk-halo decomposition 

In this Section, we decompose the total kinematics of 
the sample galaxies into the dynamical contributions of 
the baryonic disks and DM halos by disentangling the 
mass models of baryons from the total rotation curves. 
For this, we subtract in quadrature the rotation veloci¬ 
ties of the gas and stellar components from the total ro¬ 
tation velocities, and obtain implied rotation curves for 
the DM halos as shown in the figures in the Appendix 
(e.g., see the left-lower panel of Fig. A.3). For a quan¬ 
titative analysis of the DM distributi on in the gal axies, 
we fit GDM (Navarro, Frenk & White 1996b| 1997 here¬ 
after NF W) and spherica l pseu do-isothermal halo mod¬ 


els (e.g., [Begeman et al.j |1991| ), the two representative 
cusp- ana core-like halo rnodels, to the DM halo rotation 
curve, respectively. An Einasto profile which is a Sersic 
function in the context of GDM halos has been found to 
provide an equal or better fit to the halos in pure DM 
simulat ions compared to an NFW profile [Navarro et al.j 
(2004b). In this work, we use NFW profiles to quan- 
tify the DM halos of the sample galaxies for consistency 
with the previous DM mas s modelling of THINGS and 
simulated dwarf galaxies in Oh et al. (2011a 2011b). A 
relative comparison of the two halo models in terms of fit 
quality enables us to examine which model best describes 
the DM component of the sample galaxies, especially to¬ 
wards the centers of the galaxies. 


3.4.1. Cusp-like halo model 

The cosmologically motivated NFW halo model, the 
so-called ‘universal density profile’ which describes the 
cusp-like radial DM distribution found in DM-only 
AGDM simulations is given as, 

Pi 


Pnfw( 7 ?) = 


(i?/i?«)(l + i?/i?,)2’ 


( 6 ) 


where pi is correlated with the mean density of the Uni¬ 
verse at the time of the collapse of the halo and Rg is 
the ch aracteristic radius of the DM halo (Navarro et al. 
1996bI. This profile has been widely adopted to account 
for the DM distribution which steeply increases towards 
the centers of the halos in the simulations. This, a so- 
called cusp feature, can be well approximated by a power 
law, with 1.0 near the central region of the 

halos, giving a DM halo rotation velocity as follows. 


Vnfw( 7 ?) — V2001 


/ln(l -|- cx) — cx/ (1 -I- cx) 
a;[in(l c) - c/(l + cj] 


(7) 


where c is the concentration parameter defined as 
7 ?20 o/7?s- V 200 is the rotation velocity at a radius i? 2 oo 
where the mass density contrast with the critical density 
of the Universe exceed s 200, and x is defined as R/R200 
(Navarro et al. 1996b). In particular, the concentration 
parameter c is useful for quantifying the degree of DM 
concentration in galaxies. The measurement of the value 
of c in nearby galaxies provides an observational con¬ 
straint o n the central cusps pred icted fro m AGDM sim¬ 
ulatio ns ([McGaugh et al. 2007 see also de Blok et al. 
^003|. 


3.4.2. Core-like halo model 

As an alternative way to describe the DM distribution 
in a galaxy observationally motivated, spherical pseudo- 
isothermal halo models with a central constant-density 
core have been used in studies of galaxy rotation curves. 
The form of this core-like halo model is given as follows: 

Po 


PisoiR) = 


1 + {R/Rcy^ 


( 8 ) 


where pg and Rc are the core-density and core-radius 
of a halo, respectively. This halo model is employed to 
describe the mass distribution of a DM halo with a size¬ 
able constant density-core (pcxpo)- Similarly, the corre¬ 
sponding rotation velocity to the pseudo-isothermal halo 
potential is given by, 


^iSo(7?) = 


\ 


AirOpgR^ 


1-^atanfA' 
R \Rc ^ 


(9) 


In order to quantify the DM distribution in the sample 
galaxies, we fit the model rotation velocities of the two 
halo models as given in Eqs. and to the kinematic 
DM signature of the galaxies. As shown in the figures 
in the Appendix (e.g., the left-lower panel of Fig. A.3), 
the core-like pseudo-isothermal halo models are mostly 
preferred over the cusp-like NFW models in terms of the 
quality of the fit (i.e., based on reduced values). The 
fitted parameters of the halo models are presented in 
Table 12 

As given in Table [2 we derive the dynamical masses 
(Mdyn) of the sample galaxies using Uiso(7?max) and 
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Rmux at which the outermost part of the rotation curve 
is measured. Given that rotation curves at the adopted 
Rmux are mostly rising in our sample galaxies, the de¬ 
rived dynamical mass, M^yn with V[so(.Rmax) and i?max 
gives a lower limit of the halo mass. For reference, we 
also estimate the halo masses, M 200 of our sample galax¬ 
ies assuming an NFW halo model. However, the fitted 
values of c are unphysical (i.e., negative or close to zero) 
for a large fraction (11/26 galaxies) of the sample galax¬ 
ies. We therefore derive a c value using an empirical 
relationshi p between c — V 2 0 0 from the WMAP^® obser¬ 
vations in McGaugh et al. (2007) by substituting V 200 


with the Viso(l^max) adopts We then fit the NFW 
halo model after fixing c and leaving only V 200 as a free 
parameter, assuming a minimum disk (where the rota¬ 
tion curve is attributed to the DM halo only and the 
dynamical contribution of baryons is ignored). These c 
and V 200 values are given in brackets in Table Lastly, 
we derive the resulting halo ma ss, M 2 nn w i th the newly 
estimated V 200 using Eq. 3 in Oh et al. (|2011a[). As 


presented in Table M 200 values of the sample galax¬ 
ies are larger than M^yn values, which implies that our 
observations most likely do not reach the flat part of the 
rotation curves. This is consistent with the fact that the 
rotation curves of most sample galaxies are still rising at 
the last measured points. 

Except in those few cases, such as DDO 70, DDO 101, 
DDO 154, DDO 210, DDO 216, and Haro 36 where CDM 
NFW halo models provide comparable fits to the DM 
rotation curves, the fitted values of the NFW halo pa¬ 
rameters are unphysical (i.e., negative concentration pa¬ 
rameter c, unphysically large values of V 20 o)- This is 


., de Blok & Bosma 

2002 

|van Eymeren et al.| 



Oh et al. 2011bI. The slowly increasing DM rotation 


curves m the inner region of the sample galaxies reflect a 
halo whose gravitational potential is not deep enough to 
sustain the power-law DM density cusps that are as steep 
as poci?“^ °. Hence, the cosmologically motivated cusp¬ 
like halo models are not able to adequately describe the 
observed solid-body rotation curves of our sample dwarf 
galaxies. 

It may be argued that those galaxies that are equally 
well fitted by the two halo models indicate the possi¬ 
bility of a kinematic signature of central cusps in dwarf 
galaxies being consistent with ACDM simulations. De¬ 
spite the high-resolution of LITTLE THINGS VLA Hi 
observations, it is, however, most likely that the cen¬ 
tral regions of the galaxies are not fully sampled with 
a sufficient number of independent synthesized beams 
needed for distinguishing the inner steepness of the two 
halo models, accurate v. As quantified in the parameter, 
in Table ^ this is mainly due to the small 
size of the rotating disx (e.g., DDO 210) or the relatively 
large distance (e.g., Haro 36). Given that the difference 
between the cusp- and core-like halo models is the most 
prominent in the central regions of galaxies, higher sam¬ 
pling of the inner regions of the galaxies is required before 
making a firmer conclusion on the signature of the po¬ 
tential central cusps. We will discuss this matter in a 


more quantitative way in the following Section. 

4. DARK MATTER DISTRIBUTION 

In this Section, we compare the inner shape of ro¬ 
tation curves and DM density profiles of the sample 
galaxies with those of simulated dwarf galaxies from N- 
body-bSPH ACDM simulations in order to examine their 
DM distributions near the centers. These comparisons 
between observations and simulations allow us to esti¬ 
mate the degree of cuspiness of the central DM distribu¬ 
tion in a qualitative way. In addition, we also measure 
the inner DM density slopes of the galaxies to make a 
more quantitative comparison to the simulations. 

4.1. Rotation curve shape 

The characteristic shape of the steeply rising rotation 
curve inherent in the cusp-like DM distribution near the 
centers of simulated dwarf galaxies based on the ACDM 
paradigm can be used for a qualita tive test of the sim¬ 


ulations (iHay ashi fc Navarro 2006 see also Oh et al. 
2011b 2011a|). This qualitative test is particularly use¬ 


ful in that a direct comparison between the observed and 
predicted DM rotation curves can be made without any 
additional assumption on the shape of the DM halo (e.g., 
a spherical or triaxial halo potential) which is needed for 
converting rotation curves to the corresponding density 
profiles, and the associated additional uncertainties this 
might introduce. 

This is done in the left panel of Figs. and where 
we scale the rotation curves of both our sample galaxies 
and ACDM NFW halos with respect to the rotation ve¬ 
locity Vb .3 at a radius i?o .3 which is where the logarithmic 
slope of the cu rve is dlogV/dlogR = 0.3 (see Hayashi &| 
Navarro 2006). This enables us not only to make a rel- 


ative comparison of the rotation curves between obser¬ 
vations and simulations but also to accentuate the inner 
rotation curve shape which is sensitive to the degree of 
central DM concentration. 

All the scaled (DM only) rotation curves of the 21 LIT¬ 
TLE THINGS sample galaxies for which Spitzer 3.6/im 
image is available are overplotted in the upper-left panel 
of Fig. We also overplot the median values of the rota¬ 
tion curves in each O.IR/R0.3 bin. In the lower-left panel 
of Fig. ^ we also overplot the scaled rotation curves of 
seven dwarf galaxies from THINGS (three of them are 
also in LITTLE THINGS) as well as the two simulated 
dwarf galaxies presented in Governato et al. (2010) which 
all show a linear increase in their inner regions to the 
median values of the LITTLE THINGS rotation curves. 
In particular, the two simulated dwarf galaxies were af¬ 
fected by baryonic feedback processes (mainly repeated 
gas outflows driven by SN explosions) in such a way that 
the central c usps predicted from D M-only simulations 
are flattened (Governato et al. 2010). The flattened DM 
distribution results in slowly increasing rotation curves 


in the in ner region of the simulated galaxies (see Oh et al. 
(2011a) for further discussionV 


The Wilkinson Micr owave Anisotropy Probe (|Spergel et al.| 
|2003[ [Spergel et al.||2007| 


In line with the results in Oh et al. (2011a), the in¬ 
ner shape of the scaled rotation curves ot the LITTLE 
THINGS sample galaxies falls mostly below that of the 
DM-only NFW halo models, indicating a shallower DM 
distribution near the centers. As shown in the lower- 
left panel of Figs. Hi and 1^ the discrepancy with NFW 
models becomes more pronounced when comparing these 
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Fig. 2.— Left panels: The rotation curves of the first nine galaxies of the 26 LITTLE THINGS (in dynamical mass order) which 
are all scale d with respect to the rotati on velocity Vb.3 3-t ^0.3 where the logarithmic slope of the rotation curve is d\ogV/d\ogR = 0.3 as 
described in|Hayashi & Navarro](|2006f. The upper (DM + baryons) and lower (DM only) panels show the ones including and excluding 
the dynamical contribution by baryons, respectively. The grey solid and black solid lines with small dots indicate the CDM NEW dark 
matter rotation curves with V 200 which is > 90 kms“^ and < 90 kms“^, respectively. The dashed lines (denoted as ISO) show the best 
fitted pseudo-isothermal halo models to the galaxies. Right panels: The corresponding dark matter density profiles derived using the 
scaled rotation curves in the left panels. The grey (V 200 > 90 kms“^) and black solid lines with small dots (T 200 < 90 kms“^) represent 
the CDM NEW models with the inner density slope Q!~—1.0. The dashed lines indicate the best fitted pseudo-isothermal halo models with 
a~0.0. See Section |4] for more details. 


models to DM rotation curves derived after subtracting 
the contribution from baryons from the total kinematics 
for each of the sample galaxies. Instead, they are more in 
line with those of both the THINGS dwarf galaxies and 
simulations (DGl and DG2) where the effect of bary- 
onic feedback processes is included. The linearly (or less 
steeply) rising rotation curves in the inner region of the 
galaxies indicate a nearly constant or shallower mass dis¬ 
tribution towards the centers as found in the majority of 
nearby dwarf galaxies. The LITTLE THINGS sample 


galaxies give no clear indication of the central cusps in 
their DM halos unlike the predictions from AGDM DM- 
only simulations. 

4.2. Dark matter density profiles 

As a more direct way to examine the central DM dis¬ 
tribution of the sample galaxies and compare them with 
AGDM simulations, we derive their DM density profiles 
from the DM rotation curves decomposed in Section 
A direct conversion of the rotation curve to the corre- 
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Fig. 3. — Left panels: The rotation curves of the other nine galaxies of the 26 LITTLE THINGS (in dynamical mass order) which 
are all scale d with respect to the rotati on velocity Vb.3 ^0.3 where the logarithmic slope of the rotation curve is dlogV/dlogR = 0.3 as 
described in|Hayashi & Navarro](|2006||. The upper (DM + baryons) and lower (DM only) panels show the ones including and excluding 
the dynamical contribution by baryons, respectively. The grey solid and black solid lines with small dots indicate the CDM NFW dark 
matter rotation curves with V 200 which is > 90 kms~^ and < 90 kms~^, respectively. The dashed lines (denoted as ISO) show the best 
fitted pseudo-isothermal halo models to the galaxies. Right panels: The corresponding dark matter density profiles derived using the 
scaled rotation curves in the left panels. The grey (V 200 > 90 kms“^) and black solid lines with small dots (V 200 < 90 kms“^) represent 
the CDM NFW models with the inner density slope q;~— 1.0. The dashed lines indicate the best fitted pseudo-isothermal halo models with 
Q~0.0. See Section]^ for more details (continued). 


spending DM d ensity profile can be made by the follow¬ 
ing formula (see de Blok et al. 2001 for more details), 


p{R) = 


1 


47rG 


y dV 



RdR 


( 10 ) 


where V is a rotation velocity observed at a radius R, 
and G is the gravitational constant. For this conversion, 
we assume a spherical halo potential which is va lid for 
most nearby galaxies (Trachternach et al. 2008). This 


method has been used for deriving DM density profiles of 
dwarf and LSB disk galaxies, and proved to be reliable 
as found in the comparison of the derived DM density 
profiles of simulat ed dwarf galaxies with their input ones 
(Oh et al. 2011a). 

In the right panel of Figs. and we present the 
derived (dark) matter density profiles of the LITTLE 
THINGS sample galaxies. In addition, we calculate the 
median values of the DM density profiles as shown in 
the upper-right panel of Fig. We also overplot the 
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Fig. 4.— Left panels: The rotation curves of the remaining eight galaxies of the 26 LITTLE THINGS (in dynamical mass order) which 
are all scale d with respect to the rotati on velocity Vq.s at Rq.z where the logarithmic slope of the rotation curve is dlogV/dlogR = 0.3 as 
described in|Hayashi Sz Navarro](|2006||. The upper (DM + baryons) and lower (DM only) panels show the ones including and excluding 
the dynamical contribution by baryons, respectively. The grey solid and black solid lines with small dots indicate the CDM NFW dark 
matter rotation curves with V 200 which is > 90 kms~^ and < 90 kms~^, respectively. The dashed lines (denoted as ISO) show the best 
fitted pseudo-isothermal halo models to the galaxies. Right panels: The corresponding dark matter density profiles derived using the 
scaled rotation curves in the left panels. The grey (V 200 > 90 kms“^) and black solid lines with small dots (V 200 < 90 kms“^) represent 
the CDM NFW models with the inner density slope q;~— 1.0. The dashed lines indicate the best fitted pseudo-isothermal halo models with 
Q!~0.0. See Section]^ for more details (continued). 


DM density profiles of the THINGS and simulated dwarf 
galaxies in the lower-right panel of Fig. The radial fluc¬ 
tuation shown in some of the DM density profiles (e.g., 
IG 1613 etc.) is largely due to fluctuations in the derived 
rotation curves at the relevant radius. These are mainly 
because of either the effect of non-circular motions, the 
noise in velocity profiles with low S/N values, or both in 
the region. 

As already implied by the solid-body like inner rotation 
curve shape of the LITTLE THINGS sample galaxies in 


Section [4T| their central DM density profiles are system¬ 
atically shallower than those of the cusp-like DM density 
profiles predicted from DM-only ACDM simulations. As 
shown in the lower-right panel of Figs. andthe dif¬ 
ference is even more pronounced in the comparison of the 
DM density profiles corrected for the baryons although 
the dynamical contribution by baryons is rather insignif¬ 
icant. The sample dirr galaxies are more consistent with 
the core-like DM density profiles (dot-dashed lines in the 
right-hand frames of Figs. and than NFW-type 
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Fig. 5.— Upper-left panel; The (DM only) rotation curves (small dots) of the 21 LITTLE THINGS (including 3 THINGS galaxies) 
for which Spitzer 3.6/im image is available. These are all scal ed with respect to the rota tion velocity Vb.3 -^0.3 where the logarithmic 
slope of the rotation curve is dlogV/ dlogR = 0.3 as described in|Hayashi & Navarro] (|2006||. The ‘x’ symbol represents the median values of 
the rotation curves in each O.IR/Rq.s bin. The error bars show the icr scatter. Ijower-left panel: The scaled rotation curves of the seven 
THINGS, and the two simulated dwarf galaxies (DGl and DG2 in|Governato et al.||2010ll which are overplotted to the median values of 
the LITTLE THINGS rotation curves. The grey solid and black solid lines with small dots indicate the CDM NEW dark matter rotation 
curves with V 200 which is > 90 kms“^ and < 90 kms“^, respectively. Right panels: The corresponding dark matter density profiles 
derived using the scaled rotation curves in the left panels. The grey {V 200 > 90 kms~^) and black solid lines with small dots (V 200 < 90 
kms“^) represent the CDM NEW models with the inner density slope —1.0. See Section[^for more details. 


profiles. This is much like the THINGS dwarf galaxies, 
and the simulated dwarfs (DGl and DG2) with baryonic 
feedback processes as shown in the lower-right panel of 
Fig.0 

We also measure the inner density slopes a of the DM 
density proHles to quantify the cuspiness of the central 
DM distribution. This yields a more quantitative com¬ 
parison between the observations and simulations. As 
shown in the figures in the Appendix (e.g., panel (f) of 


Fig. A.3), we perform a least squares fit (dotted lines) to 
the inner data p oints (grey dots) within a ‘b reak radi us’. 


As described in de Blok & Bosma (2002| see also 


et al. 2011b I, we determine a break radius of a DM den- 


Oh 


sity profile where the slope changes most rapidly in the 
inner region of the profile. Following |de Blok fc Bosnia] 
(20021, we adopt the mean difference between the slopes 
which are measured including the first data point out¬ 
side the break-radius and excluding the data point at the 
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break radius, respectively, as an errorbar Aa of the inner 
density slope. We measure the inner density slopes a, of 
the galaxies from their total matter (including both DM 
halo and baryons) as well as DM-only density profiles. 
The former, a so-called ‘minimum disk assumption’ that 
attributes the total rotation curve to the DM component 
only, gives a steeper inner density slope. Meanwhile, the 
latter where the dynamical contribution by baryons is 
subtracted from the total rotation curve allows us to ex¬ 
amine the effect of the model ® on the measured inner 
density slope. The mean values of the slopes of the 26 
LITTLE THINGS dwarf galaxies are ami„=-0.42 ± 0.21 
and pdf=—0.32 ± 0.24^® for the minimum disk and 
the model disk assumptions, respectively. 

As expected, the slopes measured assuming the mini¬ 
mum disk are slightly steeper than those measured using 
the model Tj ® disk for most sample galaxies. However, 
the difference between the two slopes is largely insignifi¬ 
cant since most of the sample galaxies are DM dominated 
as indicated by their low baryonic fraction. The mea¬ 
sured logarithmic inner density slopes a of our sample 
galaxies are listed in Table 

The LITTLE THINGS sample dwarf galaxies do not 
in general agree with the steep logarithmic inner slope 
(^1.0) of the DM density profiles predicted from ACDM 
DM-only simulations. Instead, they show a range of shal¬ 
lower slopes being consistent with a core-like DM distri¬ 
bution at the centers, which supports the previous results 
found in nearby dwarf and LSB disk galaxies. Our results 
(Q!min=—0.42 ± 0.21; aT*.pdf=—0.32 ± 0.24) are consis¬ 
tent with the mean logarithmic slope, amin=—0.29±0.07, 
of the seven THI NGS dwarf galax ies derived assuming a 
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minimum disk in Oh et al. 
bine the sample dwarf ga 


(2011b). Moreover, if we com- 
axies from LITTLE THINGS 


and THINGS which have a similar data quality and 
whose inner density slopes are derived in exactly the 
same way, the mean value of the slopes of the 29 dwarf 
galaxies is amin=—0.40 ±0.24. This shows good agree¬ 
ment within the error bars with aniin=~0-2±0.2 derived 


from LSB galaxies in de Blok et al. (2001) but a clear 
deviation from the 1.0 predicted from ACDM DM- 
only simulations. 

A few galaxies in our sample, such as DDO 101 and 
DDO 210 whose rotation curves are equally fitted by 
both NEW and pseudo-isothermal halo models appear 
to have relatively steeper slopes comp ared to the other 
ones. However, as discussed in Section [±41 this could be 
due to insufficient sampling of the dark matter density 
profiles in the inner region. The gradient of the loga¬ 
rithmic density slope, d\ogp/d\ogR gradually decreases 
towards the outer region of a galaxy, giving a steeper 
slope a. Therefore, the steeper a slope a is, the more 
data points in the outer regions are included when mea¬ 
suring the logarithmic slope of a DM density profile. This 
could, conceivably, be the case for galaxies with insuffi¬ 
cient spatial sampling. As discussed above, the insuffi¬ 
cient sampling mainly arises from either the smaller size 
of the Hi disk (e.g., DDO 210) or the larger distances of 
the galaxies (e.g., Haro 36). 

The dynamical contribution by the stellar component is in¬ 
cluded in the DM density profiles of five galaxies where no Spitzer 
3.6/rm image is available. 


As discussed in Governato et al. (2010), the constant- 
density cores observed near the centers of dwarf galaxies 
can be reconciled with simulations by taking the effect 
of baryonic feedback processes into account without a 
need for any explicit modification of the current ACDM 
paradigm. In particular, repeated gas outflows driven 
by SNe have been found to be efficient enough to redis¬ 
tribute the matter in dwarf galaxies, resulting in shal¬ 
lower D M density profiles as observed in nearby dwarf 
galaxies (Oh et al. 2011a). 

The investigation of the effect of SN-driven gas out¬ 
flows on the central DM distribution has been extended 


to low mass field dwarf galaxies using GASOLINE (Wad- 
sley et al. 2004[), a parallel SPH tree-code with mufti- 


2012 


stepping (Governato et ah 2012 Pontzen fc Governato 
|Di (Jintio et al. 2014 see afsi^ontzen &: Governato 


2014). In the simulations, the present-day stellar mass of 
galaxies ranges from 10®'® down to lO'^'® Mq where the 
energy transfer from repeated gas outflows to the DM 
component becomes inefficient. More specifically, the SN 
feedback in small halos where less than 0.03% of the total 
amount of baryons is converted into stars is less effective 
at removing the central cusps and turning the cusp-like 
DM density profiles into core-like ones. According to this 
scenario, it is expected that the central DM distribution 
in these systems remains cuspy, and the inner slopes of 
their DM density p rofiles are steeper than those of higher 
mass counterparts (Governato et al. 2012 Zolotov et al.| 
2012). It implies that the central cusps predicted from 
AGDM DM-only simulations should survive in low mass 
halos. This underlines the cosmological importance of 
low mass dwarf galaxies in the local Universe for testing 
the ‘cusp/core’ problem. 

As discussed earlier, the systematic uncertainties 
caused by low resolution radio observations are signifi¬ 
cantly reduced in the high-resolution Hi data from LIT¬ 
TLE THINGS, which allows us to derive more accurate 
rotation curves of the sample galaxies and thus their cen¬ 
tral DM distributions. In addition, the Spitzer IRAC 
3.6/xm data combined with model T®'® values based on 
stellar population synthesis models and galaxy colors 
provide more reliable stellar masses of the sample galax¬ 
ies. 

In Fig. 1^ we plot the inner density slopes a of the 
sample dwarf galaxies from both LITTLE THINGS and 
THINGS against their stellar masses M* on a logarith¬ 
mic scale. In addition, we also add those of a sample of 
T HINGS disk galaxies whose mass models were derived 


de Blok et al. 


(2008) in order to examine the a —M* re- 


la'tionship in the higher mass regime. Compared to dwarf 
galaxies, the central kinematics of disk galaxies are usu¬ 
ally dominated by a bulge component. It is therefore not 
trivial to perform a reliable disk-halo decomposition of 
the disk galaxies despite using the multi-wavelength data 
from THINGS whose data quality is comparable to that 
of LITTLE THINGS. For this reason, in Fig. we use 
the inner density slopes of the disk galaxies measured as¬ 
suming a minimum disk. As discussed in Section [4.2[ it 
should be noted that the minimum disk assumption at¬ 
tributes the total kinematics of a galaxy to the DM only, 
ignoring the contribution of baryonic components, giving 
a lower limit on the inner density slope (i.e., a steeper 
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Fig. 6. — The inner dark matter density slope a of the sample galaxies from LITTLE THINGS (filled circles) and THINGS (filled squares, 
THINGS dwarfs; open squares, THINGS disk galaxies) against their total stellar masses. The dashed line indicates the prediction 

from ACDM dark matter only simulations. The open stars and solid line represent the a vs . of the resolved halos from the ACDM 

SPH simulations with baryonic feedback processes which are measured at 500 pc and z=0 (|Governato et al.||2012| |Di Cintio et al.||2014|l. 
See Section 13 for more discussions. 


slope). 

As shown in Fig.|^ the ACDM SPH simulations includ¬ 
ing the effect of baryonic feedback processes predict cusp¬ 
like DM distributions with steeper inner density slopes 
(a < —1.0) in DM halos whose stellar mass is less than 
about 10® Mq. As discussed in Governato et al. (2012), 
the lower the mass of the stellar component m a galaxy, 
the less the dynamical effect of SN feedback on the DM 
potential. Consequently, this results in that the central 
DM distribution in the low mass halo regime remains 
cuspy. According to this, the initial cusps formed in the 
early Universe would still exist today in low mass dwarf 
galaxies which have stellar masses less than ~ 10® Mq 


where the repeated thermal energy injection from SN ex¬ 
plosions becomes substantially inefficient, mainly due to 
the rapid ly decreasing star formation efficiency in these 
syste ms (Governato et al. 2012 Pontzen & Governato 


20121 . 


However, in Fig. as the stellar mass of a galaxy in¬ 
creases, the inner DM density slope a becomes shallower 
in the ACDM SPH simulations including baryonic feed¬ 
back processes (open stars) with respect to the coun¬ 
terpart in t he DM-only simulation s (d ashed line). As 
discuss e d in|G overnato et ^(|2010|) and|Pontz en & Gov- 
ernatol (2012p (see also |Di Cintio et al.||2014|), the cen- 
cusps can be disrupted by the rapid gas injec- 
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tion caused by SN-driven gas outflows into the central 
region of galaxies, resulting in a shallower DM density 
distribution. This shows that DM-baryon interactions in 
dwarf galaxies through gas outflows play a critical role 
not only in forming bulgeless dwarf galaxies but also in 
flattening central cusps predicted from DM-only ACDM 
simulations. This demonstrates that proper modeling 
of DM-baryon interactions in hydrodynamical ACDM 
galaxy simulations is able to alleviate the long-standing 
tension associated with the central DM distribution in 


disk (~60" diameter) and distinguish between cusp- and 
core-like DM behaviour near the center of DDO 210. We 
note that DDO 210 has the most compact Hi disk (in 
terms of the beam size, e.g., -RmaxHiJCjg^jjj^O.O as given 
in Table of the sample galaxies. Likewise, some of 
our sample galaxies may still suffer from beam smearing 
which is discussed in the following section. 

6. INNER DENSITY SLOPE VS. RESOLUTION 

As m I 
tion of 


de Blok et al. (20011, for a quantitative examina- 


refer to Governato et al. (2012|) (see alsolPontzen & Gov- 

we plot the inner density slopes a of the galaxies includ- 

ernato 2012 

21)14 and Di Gintio et ai. 

2U14I) tor a detailed 

ing the THINGS sample as well as the two simulated 


distribution in dwarf galaxies. 

Meanwhile, in Fig. the trend of slope change 
dloga/dlogM* predicted from dwarf galaxy simulations 
with baryonic feedback processes is reversed in massive 
disk galaxies (open boxes) where a bulge component be¬ 
comes dominant in the central kinematics. As discussed 


Di Cintio et al. (2014), the effect of SN feedback can 


be surpassed by the deep gravitational potential which is 
caused by the bulge component in the central region. 
However, as discussed earlier, the steep slopes of the 
THINGS disk galaxies could be partially affected by the 
minimum disk assumption used for deriving their DM 
density profiles. The central kinematics of a disk galaxy 
with a substantial bulge component is sensitive to even 
small uncertainties in T* when converting the luminosity 
profile to the mass density profile. This makes it diffi¬ 
cult to perform an accurate disk-halo decomposition of 
bulge-dominated disk galaxies. This stresses the useful¬ 
ness of bulge-less dwarf galaxies in testing the effect of 
SN feedback on the central DM distribution of galaxies. 

As shown in Fig.[^ the simulated dwarf galaxies which 
have comparable stellar masses ranging from 10® to 10 ® 
Mq show good agreement with the majority of the sam¬ 
ple dwarf galaxies from THINGS and LITTLE THINGS. 
However, as already discussed in Section ^ the rota¬ 
tion curves of some LITTLE THINGS sample galaxies 
(e.g., DDO 70, DDO 101, DDO 154, DDO 210 and Haro 
36) are equally well fitted by CDM NEW and pseudo- 
isothermal halo models in terms of values. Moreover 
some of them, such as DDO 101, DDO 210 and Haro 36, 
appear to have relatively steep inner density slopes (al¬ 
though Haro 36 is defined as a blue compact dwarf galaxy 
which usually shows a steep increase in rotation velocity 
in the inner region mainly due to young starburst compo- 
nents forme d during galaxy interaction or mergers; e.g., 
Bekki 2008). In particular, DDO 210 is a good candi- 
date tor testing the SN feedback efficiency scenario by 
Governato et al. (2012) given that its stellar mass falls 
within the regime where according to the simulations pri¬ 
mordial CDM cusps are expected to survive. In Fig. 
DDO 210 shows no distinct signature of the central cusp 
(a~—1.0) given its corresponding stellar mass, deviat¬ 
ing from the prediction of ACDM DM-only simulations 
(dashed line). Nevertheless, in a qualitative sense, its 
relatively steep inner density slope 0.70 compared 
to the others still leaves a room for a potential signature 
of the central cusp, which supports the lower SN energy 
injection scenario in low mass dwarf galaxies. However, 
as noted earlier, the resolution of LITTLE THINGS Hi 
observations is not high enough to resolve the small Hi 


et al. (2010) against the observed radii of their innermost 
point /tinner in Fig.[^ For the sample dwarf galaxies from 
LITTLE THINGS, we use the slopes derived assuming 
the model T®-® disk. We also show the a-Rinner rela¬ 
tions of the NEW and pseudo-isothermal halo models as 
solid and dotted lines, respectively, derived using their 
analytical formulas as given in Eqs. |^and[^ 

As shown in Fig.[^ most sample galaxies show signifi¬ 
cant deviations from the predicted a-/tinner trend (solid 
line) of ACDM NEW halos at around a /tinner of ~0.2 
kpc. Instead, they are more consistent with those of 
pseudo-isothermal halo models with different core-radii 
(dotted lines) as well as the earlier results found from 
LSB galaxies (grey symbols) inlde Blok & Bosma (2002). 
However, the clear difference between the two halo mod- 
els (i.e., NEW and pseudo-isothermal) at high resolu¬ 
tions (e.g., /tinner<0.5 kpc) becomes ambiguous as the 
innermost radius /tinner of a given DM density profile in¬ 
creases. For example, a galaxy with a larger /tinner (he., 
low resolution) tends to show a steeper inner slope of DM 
density profile. The larger /tinner makes it lie in the re¬ 
gion where the slopes of the two halo models are approxi¬ 
mately simi lar to each other. In addition, as discussed in 
Section 4.2 the derivative dlogp/dlog/t of a DM density 
profile on a logarithmic scale decreases towards the outer 
region of a galaxy. If the DM density profile is affected 
by beam smearing, the break radius of the profile which 
is determined when measuring the inner slope tends to 
migrate into the outer regime where dlogp/dlogR has a 
lower value. Therefore, the inner density slope a within 
the break radius is most likely to be steeper than the 
ones derived from well sampled profiles. This could be 
the case of DDO 101 and DDO 210. In particular, DDO 
101 is most likely to be affected by the beam smearing 
effect as shown in Fig. 

Yet higher resolution velocity fields obtained with ra¬ 
dio interferometers or using other tracers such as integral 
field mapping are required to study the effect of SN feed¬ 
back on the central cusps of the lowest mass dwarf galax¬ 
ies. Such high-resolution observations of low mass dwarf 
galaxies would provide an ultimate test of the ACDM 
paradigm. Unlike clusters of galaxies where the depth of 
the gravitational potential well is deep enough to retain 
warm DM (WDM) as well as CDM, there is no room for 
WDM in dwarf galaxies inhabiting DM halos with much 
shallower potential wells. Therefore, finding a signature 
of a central cusp in dwarf galaxies will prove that there 
is at least some cold DM in the Universe. This again 
highlights the cosmological importance of low mass dwarf 
galaxies, not only for resolving the ‘cusp/core’ contro- 
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Fig. 7. — The inner slope of t he dark matter densi ty profiles a vs. the radius Rin of the innermost point within which a is measured 
as described in the small figure l |de Blok et al.||2001[|. The Q;-Rj„ of the sample galaxies from LITTLE THINGS, THINGS and the two 
simulated dw arf galaxies (DGl a nd DG2: i Governato et al. 1 1201011 as well as t he previous measure ments (grey symbols) of LSB galaxies 
(open circles: |de Blok et al.||2001| triangles:T3^^E^^^BosmapDDz| open stars: jSwaters et al.||2003||. Filled circles with arrows indicate the 
galaxies of which inner density slopes are measured assuming a ‘minimum disk’, giving a steeper slope. The solid and dotted lines represent 
the a-Rin trends of dark-matter-only AGDM NEW and pseudo-isothermal halo models, respectively. See Section]^ for more details. 


versy in ACDM simulations but also as an indirect proof 
for the existence of CDM in the Universe. 

7. CONCLUSION 

In this paper we derive the rotation curves of 26 dwarf 
galaxies culled from LITTLE THINGS, and examine 
their DM distributions near the centers of the galaxies. 
From this, we address the ‘cusp/core’ problem which has 
been one of the long-standing problems in ACDM simu¬ 
lations on galactic scales. The high-resolution LITTLE 
THINGS Hi data (^6"angular; ^2.6 kms“^ spectral) 
complemented with optical and Spitzer IRAC 3.6/rm im¬ 
ages are sufficiently detailed to resolve the central region 
of the sample galaxies where the cusp- and core-like halo 


models are clearly distinguished. 

In particular, we use the bulk velocity fields of the 
galax i es ex tracted using the method described in |Oh| 


et al. (2008) to correct for turbulent random non-circular 


gas motions. This enables us to derive more reliable rota¬ 
tion curves and thus more accurate DM distributions in 
the galaxies. We corrected for the modest dynamical con¬ 
tribution by baryons in dwarf galaxies by using Spitzer 
IRAC 3.6ixm images combined with model Tj ® values 
based on stellar population synthesis models. This al¬ 
lowed us to derive robust mass models of the stellar com¬ 
ponents of the galaxies and thus better constrain their 
central DM distributions. 

From this, we found that the decomposed DM rotation 
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curves of most sample galaxies are well matched in shape 
to those of core-like halos which are characterised by a 
linear increase of rotation velocity in the inner region. 
We also derive the DM density profiles of the sample 
galaxies and quantify the degree of the central DM con¬ 
centration by measuring the logarithmic inner slopes of 
the profiles. The mean value of the inner slopes a of 
the 26 sample galaxies is —0.32, which indicates a mass 
distribution with a sizable constant density-core towards 
the centers of the galaxies. This is consistent with that 
found in most nearby dwarf g alaxies (e.g., LSB ga lax- 
ies in lde Blok fc Bosma|[200^ THINGS dwarfs in |Oh| 
et al. |2011b I which all show a linear increase in the in- 
ner shapes of their rotation curves resulting in shallower 
inner density slopes (< a —0.2). Considering the 
fact that observational uncertainties are significantly re¬ 
duced in the high-resolution LITTLE THINGS data, the 
core-like DM distribution found in our sample galaxies 
provides a stringent observational constraint on the cen¬ 
tral DM distribution of halos in ACDM simulations. 

We find that the derived slopes of the DM density pro¬ 
files do not agree with the cusp predicted by ACDM DM- 
only simulations. However, recent cosmologi c al N- body 
SP H galaxy simulations by Governato et al. | (|2010|) (see 


of the galaxies acts to results in steeper observed inner 
DM density slopes. 

According to the latest N-body SPH simulations of 


__ dMg 

also I Governato et al.| |2012| |Pontzen & Governato||2012 1 

have shown that the discrepancy between observations 
and simulations can be reconciled within the ACDM 
paradigm by considering the dynamical effect of bary- 
onic feedback processes on the central cusps. Accord¬ 
ing to the simulations, DM-baryon interactions in dwarf 
galaxies through gas outflows driven by SN explosions 
play a critical role not only in forming bulgeless dwarf 
galaxies but also i n turning c e ntral c usps into cores. 

As discussed in Oh et al. (2011a), the slowly rising 
rotation curves and the resulting shallower DM density 
profiles of the simulated dwarf galaxies with SN feedback 
are qualitatively similar to those of the dwarf galaxies 
from THINGS. This is also the case for our sample galax¬ 
ies from LITTLE THINGS whose DM rotation curves in 
the inner region rise too slowly to match the steep ro¬ 
tation curves of CDM halos. This shows that proper 
modelling of DM-baryon interactions in ACDM galaxy 
simulations is able to alleviate the long-standing tension 
between observations and simulations regarding the cen¬ 
tral DM distribution in dwarf galaxies. 

Notwithstanding the dominant trend of core-like DM 
distribution in the LITTLE THINGS sample galaxies, 
some of the sample galaxies, such as DDO 210 and Haro 
36 are equally well fitted by core- and cusp-like halo mod¬ 
els in describing their DM rotation curves. They have rel¬ 
atively steeper inner density slopes with —0.70 and 
—0.50 for DDO 210 and Haro 36, respectively, compared 
to the mean value (DM only) of the rest of the sample 
(^ —0.29). It is possible, however, despite the high an¬ 
gular resolution of the LITTLE THINGS Hi data, the 
inner density slopes of these two galaxies are affected by 
beam smearing (the Hi disk of DDO 210 being intrin¬ 
sically small and Haro 36 being one of the more distant 
dwarfs). Insufficient spatial resolution in the inner region 


dwarf galax i es wit h baryonic feedback processes (Gover¬ 
nato et al. 2012), the SN feedback in low mass dwarf 


galaxies with a stellar mass less than IO^Mq is not suf¬ 
ficient to disrupt the central cusps, the repeated energy 
injection from SN explosions into DM halos becoming 
inefficient, largely due to low star formation efficiencies 
in these low mass systems. However, previous obser¬ 
vational studies regarding the ‘cusp/core’ problem have 
mostly focused on relatively massive dwarf galaxies for 
which reliable rotation curves are available. Low mass 
dwarf galaxies have been usually excluded for the study 
of the central DM distribution due to the low amplitude 
of their maximum rotation velocities. The kinematics of 
such low mass dwarf galaxies is more vulnerable to kine¬ 
matic disturbances like non-circular motions in galaxies 
compared to more massive ones, and also more sensitive 
to additional corrections (e.g., asymmetric drift) made 
when deriving rotation curves. In this respect, the pos¬ 
sibility of a selection effect in favor of relatively massive 
dwarf galaxies where the effect of SN feedback is enough 
to turn the central cusps into cores should be considered. 

It would therefore be worthwhile to perform high- 
resolution follow-up observations, for example using an 
optical integral field unit, of low mass dwarf galaxies in¬ 
cluding some of the LITTLE THINGS dwarf galaxies 
(e.g., DDO 210) whose stellar masses lie in the regime 
where primordial CDM cusps are predicted from the sim¬ 


ulations (Governato et al. 


2012). These high-resolution 


observations will enable us to achieve a finer sampling of 
the central region of the galaxies, and thus more accurate 
inner DM density profiles. From this, more stringent ob¬ 
servational constraints on the central cusp of low mass 
dwarf galaxies could be provided. Moreover, an accurate 
measurement of the DM distribution in these low mass 
dwarf galaxies will provide an ultimate test for the CDM 
paradigm given that dwarf galaxies inhabiting DM halos 
with a shallow potential well have only room for CDM, 
unlike clusters of galaxies whose gravitational potential 
is deep enough to retain WDM as well as CDM. There¬ 
fore, the presence or absence of a signature of the central 
cusp in these low mass halos will provide a critical obser¬ 
vational test, either supporting or falsifying the ACDM 
paradigm. 
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Appendix. DATA AND KINEMATIC ANALYSIS 


In this Appendix, we present the data and kinematic analysis of the 26 LITTLE THINGS (Hunter et al. 2012 1 dwarf 
galaxies. For each galaxy, we show the (1) data, (2) kinematic analysis, and (3) mass modelling with descriptions: 

A. Data— (a) Integrated Hi intensity map (moment 0). The contour levels start at +3a in steps of +3cr. (b) Intensity- 
weighted mean velocity field (moment 1). (c) Velocity dispersion map (moment 2). (d) Spitzer IRAC S.diim im age 
obtained from the archives including ^Spitzer Infrared Nearby Galaxies Survey’ (SINGS; Kennicutt et al.| 2003) and 
S pitzer ‘Local Volu me Legacy’ (LVL; [Dale et al. 2009). (e) Bulk velocity field extra c ted us ing the method described 


Oh et al. (2008). (f) Velocity field of strong uon-circular motions as in Oh et al. (2008). (g) Model velocity field 


of the tilted-r mg model derive d using the bulk velocity field in the panel (e). (h) Ve^ 


_ _ _ * . , , , ocity field of weak non-circular 

Oh et al. ( |2008 ). (i), (j) Position-velocity diagram taken along the average position angle of the major 


motions as in 

and minor axes as giveri in I'able ^ The dashed lines indicate the systemic velocity and position of the kinematic 
center derived in this paper. The bulk (black dots) and asymmetric drift corrected bulk (yellow) rotation curves are 
overplotted. The curves are converted back to radial velocities using the geometrical parameters determined by the 
tilted-ring analysis as listed in Table Moment maps and velocity fields are extracted using t he robust-weighted d ata 
cubes. The beam size is indicated by the ellipse in the bottom-right corner of each panel. See Hunter et al. (2012) for 
a detailed description of the data cube. 

B. Kinematic analysis 

• Rotation curves— The tilted-ring model derived using the bulk velocity field as given in the panel (e) of A. 
Data. The open grey circles shown in all panels are the fit results with all ring parameters (i.e., XPOS, YPOS, 
VSYS, PA, INCL and VROT) free. The grey filled dots in the VROT panel indicate the final rotation velocity derived 
using the entire velocity field after fixing other ring parameters to the values as shown in other panels (solid lines). 
The upright and upside-down triangles show the rotation velocities derived using the receding and approaching sides, 
respectively, while keeping other ring parameters at the values indicated in the other panels with solid lines. 

• Asymmetric drift correction— (a) grey dots represent the radial asymmetric drift correction, ctd- The open 
circles indicate the derived rotation velocity from the tilted-ring analysis, and the black dots show the one corrected 
for asymmetric drift, (b) Azimuthally averaged Hi velocity dispersion, (c) Radial Hi surface density derived applying 
the kinematic geometry from t he tilte d-ring analysis, (d) The dashed line shows the fit of the analytical function given 
in Eq. I^to Ehict^. See Section 3.1.2 for more details. 

• Harmonic analysis— Harmonic decompositions of the bulk (black dots) and intensity-weighted mean (moment 
1) velocity fields. The decomposition is made using the RESWRI task in GIPSY. The so lid and dashed lines in the 
midd le-right panel indicate average global elongations of the halo potential as described in Schoenmakers et al.| (1997 
19991 derived using the bulk and moment 1 velocity fields, respectively. The median absolute amplitudes ((A 
the phases {(pi, (p 2 , and ^ 3 ) of each component are presented in the bottom panels. 

C. Mass modelling 

• Mass models of baryons— (a) Spitzer IRAC 3.6/im surface brightness profile derived applying the kinematic 

geometry from the tilted-ring analysis, (b) T he stellar mass-to-light value in the 3.6/im, ® derive d using the 


empirical relat ion described in Oh et al. 

Charlotl (|2003|) andlHell & de Jongl (|200ip. We refer to 

t1i 


(2008) which is based on ste llar population synthesis models in Bruzual & 
Oh et al.| (|2008|) for a full description, (c) Stellar miass surface 
11 


density prohle in the 3.6/im. (d) ine rotation velocity for the stellar component derived using the mass density profile 
in panel (c). (e) Azimuthally averaged optical color, B — V. (f) Gas column density derived applying the kinematic 
geometry from the tilted-ring analysis, (g) Mass surface density profile for the gas component which is scaled up by 
1.4 to account for Helium and metals, (h) The resulting rotation velocity for the gas component. 

• Disk-halo decomposition— The grey dots indicate the bulk rotation curve which is corrected for asymmetric 
drift. The dotted and dash-dotted lines represent the rotation velocities of the gas and stellar components, respectively. 
The open circles show the rotation velocity of dark matter halo only, derived after subtracting the contribution to 
the rotational velocity of the baryons. The dashed and solid lines are the fits of ACDM NEW and pseudo-isothermal 
halo models to the dark matter-only curve (open circles). The reduced value for each halo model is denoted. The 
dots (DM rotation curve — NEW model) and open circles (DM rotation curve — pseudo-isothermal model) in the 
lower panel represent the velocity differences between the rotation curve of the dark matter halo and the best fit halo 
models. 

• Mass density profile— The circles indicate the mass density profile derived from the asymmetric drift corrected 
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bulk rotation curve which includes the dynamical contributions by dark matter halo and baryons (i.e., minimum disk 
assumption). The squares show the one derived from the dark matter rotation curve where the dynamical contribution 
by baryons is subtracted (dark matter only). The logarithmic inner slope a of the dark matter density profile is 
measured by a least squares fit (dotted line) to the inner data points (filled squares). The dashed and solid lines show 
the density profiles derived from the best fit NFW and pseudo-isothermal halo models, respectively. 
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Rotation curves, asymmetric drift correction and harmonic analysis of CVnIdwA. See Appendix section B for details. 
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Fig. a. 3.— The mass models of baryons, disk-halo decomposition and mass density profile of CVnIdwA. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.4.— Hi data and Spitzer IRAC S.Qpm image of DDO 43. The systemic velocity is indicated by the thick contours in the velocity 
fields, and the iso-velocity contours are spaced by 10 kms“^. Velocity dispersion contours run from 0 to 20 kms~^ with a spacing of 5 
kms~^. See Appendix section A for details. 
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Fig. a. 5.— Rotation curves, asymmetric drift correction and harmonic analysis of DDO 43. See Appendix section B for details. 
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Fig. a. 6.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 43. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.7.— Hi data and Spitzer IRAC S.Qpm image of DDO 46. The systemic velocity is indicated by the thick contours in the velocity 
fields, and the iso-velocity contours are spaced by 10 kms“^. Velocity dispersion contours run from 0 to 20 kms~^ with a spacing of 5 
kms~^. See Appendix section A for details. 

























Phase 

(radian) Velocity (kms"') Velocity (km s"') VSYS (km s“') VROT (km s" 


30 


Oh et al. 


DDO 46 Rotation Curves 






Asymmetric Drift Correction 



0123 0123 0123 0123 

R (kpc) R (kpc) R (kpc) R (kpc) 


Harmonic Analysis 






0 1 2 
R (kpc) 


3 0 12 

R (kpc) 


3 0 12 3 

R (kpc) 


2 3 

m 


CM 


Fig. a.8 


Rotation curves, asymmetric drift correction and harmonic analysis of DDO 46. See Appendix section B for details. 
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Fig. a. 9.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 46. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.10.— Hi data and Spitzer IRAC image of DDO 47. The systemic velocity is indicated by the thick contours in the velocity 

fields, and the iso-velocity contours are spaced by 10 kms“^. Velocity dispersion contours run from 0 to 20 kms~^ with a spacing of 5 
kms~^. See Appendix section A for details. 
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Fig. a. 11.— Rotation curves, asymmetric drift correction and harmonic analysis of DDO 47. See Appendix section B for details. 
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Fig. a. 12.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 47. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.13.— Hi data and Spitzer IRAC image of DDO 50. The systemic velocity is indicated by the thick contours in the velocity 

fields, and the iso-velocity contours are spaced by 15 kms“^. Velocity dispersion contours run from 0 to 30 kms“^ with a spacing of 20 
kms~^. See Appendix section A for details. 
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Fig. a. 14.— Rotation curves, asymmetric drift correction and harmonic analysis of DDO 50. See Appendix section B for details. 
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Fig. a. 15.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 50. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a. 17.— Rotation curves, asymmetric drift correction and harmonic analysis of DDO 52. See Appendix section B for details. 
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Fig. a. 18.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 52. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a. 20.— Rotation curves, asymmetric drift correction and harmonic analysis of DDO 53. See Appendix section B for details. 
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Fig. a. 21.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 53. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.22.— Hi data and Spitzer IRAC SSp^m image of DDO 70. The systemic velocity is indicated by the thick contours in the velocity 
fields, and the iso-velocity contours are spaced by 10 kms“^. Velocity dispersion contours run from 0 to 30 kms“^ with a spacing of 10 
kms~^. See Appendix section A for details. 
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Fig. a. 23.— Rotation curves, asymmetric drift correction and harmonic analysis of DDO 70. See Appendix section B for details. 


<A> ep„,sin20, 
(km s”') 







































































Velocity (km s ') orcsec 


46 


Oh et al. 


DDO 70 


Mass models of boryons 




J M 


R (kpc) 



o 

Q. 


: M 




1 2 
R (kpc) 



R (kpc) 


Disk —holo decomposition Moss density profile 




Fig. a. 24.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 70. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.25.— Hi data and Spitzer IRAC SSp^m image of DDO 87. The systemic velocity is indicated by the thick contours in the velocity 
fields, and the iso-velocity contours are spaced by 10 kms“^. Velocity dispersion contours run from 0 to 20 kms~^ with a spacing of 5 
kms~^. See Appendix section A for details. 
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Fig. a. 26.— Rotation curves, asymmetric drift correction and harmonic analysis of DDO 87. See Appendix section B for details. 


<A> ep„,sln20, 
(km s”') 









































































































Velocity (km s ') orcsec 


High-resolution mass models of dwarf galaxies from LITTLE THINGS 


49 


DDO 87 


Mass models of boryons 



0.5 

r (b)i 

^ 6 

CN 

- (c)- 

30 

r (d)i 

0.4 


k 4 


'm 20 

L _I 

0.3 


CN 

°IA|) 

0 

“ @ ” 

£ 

^ 10 


0.2 


0 

U) 

M 0 


1 0 



0 12 3 4 5 6 7 


0 12 3 4 5 6 7 


0 12 3 4 5 6 7 



o 

Q. 




R (kpc) 



0 1 2 3 4 5 6 7 

R (kpc) 


Moss density profile 



Fig. a. 27.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 87. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.28. — Hi data and Spitzer IRAC 3.6/im image of DDO 101. The systemic velocity is indicated by the thick contours in the velocity 
fields, and the iso-velocity contours are spaced by 10 kms“^. Velocity dispersion contours run from 0 to 20 kms“^ with a spacing of 10 
kms~^. See Appendix section A for details. 
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Fig. a. 29.— Rotation curves, asymmetric drift correction and harmonic analysis of DDO 101. See Appendix section B for details. 
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Fig. a. 30.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 101. Please refer to the text in 
Sections [ 3 ] and O for full information. 
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Fig. a. 31. — Hi data and Spitzer IRAC 3.6/im image of DDO 126. The systemic velocity is indicated by the thick contours in the 
velocity fields, and the iso-velocity contours are spaced by 8 kms“^. Velocity dispersion contours run from 0 to 20 kms~^ with a spacing 
of 5 kms“^. See Appendix section A for details. 
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Fig. a. 32.— Rotation curves, asymmetric drift correction and harmonic analysis of DDO 126. See Appendix section B for details. 
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Fig. a. 33.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 126. Please refer to the text in 
Sections [ 3 ] and O for full information. 
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Fig. a. 34. — Hi data and Spitzer IRAC 3.6/im image of DDO 133. The systemic velocity is indicated by the thick contours in the 
velocity fields, and the iso-velocity contours are spaced by 10 kms~^. Velocity dispersion contours run from 0 to 20 kms“^ with a spacing 
of 5 kms“^. See Appendix section A for details. 
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Fig. a. 35.— Rotation curves, asymmetric drift correction and harmonic analysis of DDO 133. See Appendix section B for details. 
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Fig. a. 36.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 133. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.37.— Hi data and Spitzer IRAC 3.6/im image of DDO 154. The systemic velocity is indicated by the thick contours in the velocity 
fields, and the iso-velocity contours are spaced by 10 kms“^. Velocity dispersion contours run from 0 to 20 kms“^ with a spacing of 15 
kms~^. See Appendix section A for details. 
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Fig. a. 38.— Rotation curves, asymmetric drift correction and harmonic analysis of DDO 154. See Appendix section B for details. 
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Fig. a. 39.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 154. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a. 40. — Hi data and Spitzer IRAC 3.6/im image of DDO 168. The systemic velocity is indicated by the thick contours in the 
velocity fields, and the iso-velocity contours are spaced by 10 kms~^. Velocity dispersion contours run from 0 to 20 kms“^ with a spacing 
of 5 kms“^. See Appendix section A for details. 
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Fig. a. 41.— Rotation curves, asymmetric drift correction and harmonic analysis of DDO 168. See Appendix section B for details. 
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Fig. a. 42.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 168. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.43.— Hi data and Spitzer IRAC 3.6/im image of DDO 210. The systemic velocity is indicated by the thick contours in the 
velocity fields, and the iso-velocity contours are spaced by 2 kms“^. Velocity dispersion contours run from 0 to 20 kms~^ with a spacing 
of 2 kms“^. See Appendix section A for details. 
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Fig. a. 44.— Rotation curves, asymmetric drift correction and harmonic analysis of DDO 210. See Appendix section B for details. 
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Fig. a. 45.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 210. Please refer to the text in 
Sections [ 3 ] and O for full information. 
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Fig. a. 46. — Hi data and Spitzer IRAC 3.6/im image of DDO 216. The systemic velocity is indicated by the thick contours in the 
velocity fields, and the iso-velocity contours are spaced by 8 kms“^. Velocity dispersion contours run from 0 to 20 kms~^ with a spacing 
of 4 kms“^. See Appendix section A for details. 
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Fig. a. 47.— Rotation curves, asymmetric drift correction and harmonic analysis of DDO 216. See Appendix section B for details. 
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Fig. a. 48.— The mass models of baryons, disk-halo decomposition and mass density profile of DDO 216. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.49.— Hi data and Spitzer IRAC 3.6^m image of F564-V3. The systemic velocity is indicated by the thick contours in the velocity 
fields, and the iso-velocity contours are spaced by 8 kms~^. Velocity dispersion contours run from 0 to 20 kms~^ with a spacing of 10 
kms~^. See Appendix section A for details. 
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Fig. a. 50.— Rotation curves, asymmetric drift correction and harmonic analysis of F564-V3. See Appendix section B for details. 
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Fig. a. 51.— The mass models of baryons, disk-halo decomposition and mass density profile of F564-V3. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a. 53.— Rotation curves, asymmetric drift correction and harmonic analysis of IC 10. See Appendix section B for details. 


<A> ep„,sin20, 
(km s”') 









































































































































Velocity (km s ') At3.6('^*3g orcsec 


76 


Oh et al. 


1C 10 


Mass models of boryons 




^11 

o 


50 
0 


15 


^ 10 
o 



0.2 0.4 0.6 

R (kpc) 


0.2 0.4 0.6 

R (kpc) 


0.2 0.4 0.6 

R (kpc) 


Disk-holo decomposition 



Moss density profile 



R (kpc) 


Fig. a. 54.— The mass models of baryons, disk-halo decomposition and mass density profile of IC 10. Please refer to the text in Sections]^ 
and m for full information. 
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Fig. a.55.— Hi data and Spitzer IRAC 3.6/im image of IC 1613. The systemic velocity is indicated by the thick contours in the velocity 
fields, and the iso-velocity contours are spaced by 10 kms“^. Velocity dispersion contours run from 0 to 20 kms“^ with a spacing of 10 
kms~^. See Appendix section A for details. 
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Fig. a. 56.— Rotation curves, asymmetric drift correction and harmonic analysis of IC 1613. See Appendix section B for details. 
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Fig. a. 57.— The mass models of baryons, disk-halo decomposition and mass density profile of IC 1613. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.58.— Hi data and Spitzer IRAC 3.6/im image of NGC 1569. The systemic velocity is indicated by the thick contours in the 
velocity fields, and the iso-velocity contours are spaced by 25 kms~^. Velocity dispersion contours run from 0 to 20 kms~^ with a spacing 
of 10 kms“^. See Appendix section A for details. 
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Fig. a. 59.— Rotation curves, asymmetric drift correction and harmonic analysis of NGC 1569. See Appendix section B for details. 
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Fig. a. 60.— The mass models of baryons, disk-halo decomposition and mass density profile of NGC 1569. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.61.— Hi data and Spitzer IRAC 3.6/im image of NGC 2366. The systemic velocity is indicated by the thick contours in the 
velocity fields, and the iso-velocity contours are spaced by 20 kms~^. Velocity dispersion contours run from 0 to 20 kms~^ with a spacing 
of 20 kms“^. See Appendix section A for details. 
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Fig. a. 62.— Rotation curves, asymmetric drift correction and harmonic analysis of NGC 2366. See Appendix section B for details. 
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Fig. a. 63.— The mass models of baryons, disk-halo decomposition and mass density profile of NGC 2366. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.64.— Hi data and Spitzer IRAC 3.6/im image of NGC 3738. The systemic velocity is indicated by the thick contours in the 
velocity fields, and the iso-velocity contours are spaced by 10 kms~^. Velocity dispersion contours run from 0 to 20 kms~^ with a spacing 
of 10 kms“^. See Appendix section A for details. 
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Fig. a. 65.— Rotation curves, asymmetric drift correction and harmonic analysis of NGC 3738. See Appendix section B for details. 
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Fig. a. 66.— The mass models of baryons, disk-halo decomposition and mass density profile of NGC 3738. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.67.— Hi data and Spitzer IRAC 3.6/im image of UGC 8508. The systemic velocity is indicated by the thick contours in the 
velocity fields, and the iso-velocity contours are spaced by 10 kms~^. Velocity dispersion contours run from 0 to 20 kms“^ with a spacing 
of 10 kms“^. See Appendix section A for details. 
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Fig. a.68 .— Rotation curves, asymmetric drift correction and harmonic analysis of UGC 8508. See Appendix section B for details. 
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Fig. a. 69.— The mass models of baryons, disk-halo decomposition and mass density profile of UGC 8508. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 
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Fig. a.70.— Hi data and Spitzer IRAC 3.6/im image of WLM. The systemic velocity is indicated by the thick contours in the velocity 
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Fig. a. 71.— Rotation curves, asymmetric drift correction and harmonic analysis of WLM. See Appendix section B for details. 
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Fig. a. 72.— The mass models of baryons, disk-halo decomposition and mass density profile of WLM. Please refer to the text in Sections]^ 
and m for full information. 
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Fig. a.73.— Hi data and Spitzer IRAC 3.Gpm image of Haro 29. The systemic velocity is indicated by the thick contours in the velocity 
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Fig. a. 74.— Rotation curves, asymmetric drift correction and harmonic analysis of Haro 29. See Appendix section B for details. 
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Fig. a. 75.— The mass models of baryons, disk-halo decomposition and mass density profile of Haro 29. Please refer to the text in 
Sections [ 3 ] and O for full information. 
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Fig. a.76.— Hi data and Spitzer IRAC 3.6/i.m image of Haro 36. The systemic velocity is indicated by the thick contours in the velocity 
fields, and the iso-velocity contours are spaced by 8 kms“^. Velocity dispersion contours run from 0 to 20 kms“^ with a spacing of 5 
kms~^. See Appendix section A for details. 
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Fig. a. 77.— Rotation curves, asymmetric drift correction and harmonic analysis of Haro 36. See Appendix section B for details. 


<A> ep„,sin20, 
(km s”') 





















































Velocity (km s ') orcsec 


100 


Oh et al. 


Haro 36 


Mass models of baryons 



0.5 

- (b)- 

^ 20 

- (c)- 

30 

- (d)- 

0.4 


1 

u 

Q. 

S 

'm 20 




10 

0 J 

£ 

; ; 

0.3 


'w' 


^ 10 

- -- 

0.2 


0 

^ r\ 

M 0 

0 

_ ®©©©o®®- 

0 

s> 

0 

L J 








0123 0123 0123 


E 

o 

(/> 

£ 

o 

o 


w 





R (kpc) 


Disk-holo decomposition 


Moss density profile 




R (kpc) 


Fig. a. 78.— The mass models of baryons, disk-halo decomposition and mass density profile of Haro 36. Please refer to the text in 
Sections [ 3 ] and [ 4 ] for full information. 



































